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Abstract 

We study the electronic and magnetic structure of carbon and nitrogen impurities and interstitials 
in rutile Ti02. To this end we perform ab-initio calculations of a 48-atom supercell employing the 
VASP code. In order to obtain a realistic description of the electronic and magnetic structure, 
exchange and correlation are treated with the HSE06 hybrid functional. Both, atomic positions 
and cell dimensions are fully relaxed. Substitutional carbon and nitrogen are found to have a 
magnetic moment of 2 and IpB, respectively, with a tendency for anti-ferromagnetic long range 
order. For C/N on interstitial sites we find that carbon is non-magnetic while nitrogen always 
possesses a magnetic moment of IpB- We find that these interstitial positions are on a saddle 
point of the total energy. The stable configuration is reached when both carbon and nitrogen form 
a C-0 and N-0 dimer with a bond length close to the double bond for CO and NO. This result 
is in agreement with earlier experimental investigations detecting such N-0 entities from XPS 
measurements. The frequencies of the symmetric stretching mode are calculated for these dimers, 
which could provide a means for experimental verification. For all conhgurations investigated both 
C and N states are found inside the Ti02 gap. These new electronic states are discussed with 
respect to tuning doped Ti 02 for the application in photocatalysis. 

PACS numbers: 61.72.uf, 61.72.uj - doping of Semiconductors 

71.15.Mb - Density functional theory, local density approximation, gradient and other 
75.50.Pp - magnetic Semiconductors 
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I. INTRODUCTION 


Ti02 plays an important role in the fields of dilute magnetic semiconductors (DMS) and 
photocatalysis. For both applications Ti02 requires tailor made electronic properties, which 
can be manipulated by selectively diluting the system by insertion of dopants. 

DMS have been investigated in the past in order to gain new insights in the functional 
principles and the manipulation of the magnetic and semiconducting properties. The search 
for magnetic semiconductors or for half-metals used for spin-injection led to the discovery 
of the new material class of p-electron magnets. The prospect of controlling the charge 
and the electron spin as information carriers would make it possible to combine information 
processing and storage at the same time.^^^ The second major application is in the field 
of photocatalysis, where the favourable electronic gap range of 1.65 — 3.1 eV can be also 
achieved by doping^^® . 

Regarding rutile Ti02 from a photocatalytic point of view we are especially interested 
in the manipulation of the band gap by incorporation of impurities. Rutile Ti02 has a 
band gap of 3.0 eV, which corresponds to an absorption spectrum in the ultraviolet region. 
The absorption of photons lead to the creation of electron hole pairs, which are of special 
relevance for photocatalytic and photoelectrochemical applications like the photo-induced 
decomposition of water on Ti02 electrodes demonstrated by Fujishima and Honda.^ Another 
line of research regards the photocatalytic properties; due do the creation of electron hole 
pairs the valence band (VB) becomes oxidative and the conduction band (CB) reductive. 
Surrounding molecules can therefore be oxidized or reduced and as a consequence form 
radicals which are harmful to organic compounds like bacteria and fungi. Hence, Ti02 can 
be used as a purifier, removing gaseous or aqueous contaminants.®’® Beside water- and air- 
purification there is a wide range of other photocatalytic applications, like antifogging and 
self-cleaning surfaces.^® Since UV light only makes up 5% of the sunlight spectrum one aims 
to shift the absorption spectrum of Ti 02 into the visible region to improve the absorption 
rate and therefore the photocatalytic efficiency. Several studies have shown that impurities 
like nitrogen and carbon verifiable modify the band gap. Beside computational studies done 
by DiValentin et al.^^’^^ and Yang et al.^® there are experimental studies confirming the 
earlier theoretical results. Diwald et al.^"^ reported a blueshift in the band gap after doping 
with nitrogen. In a second study^® they observed a redshift by inducing nitrogen into the 
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Ti02 rutile host matrix. For the first case they assumed oxygen atoms being substituted 
by nitrogen impurities and for the latter nitrogen atoms being placed on interstitial sites. 
Batzill et al.^® on the other hand reported to see a redshift when substituting O by N. 
Motivated by these partly contradictory reports we perform our theoretical investigation in 
order to interpret the experimental results. 

In the present study we investigate Ti02 rutile doped with carbon and nitrogen atoms, 
respectively, where we assume C and N either to replace oxygen or to occupy interstitial 
positions. Although magnetic order is most common in metallic materials with narrow 
bands of d- or /- electrons, the carriers of the magnetic moments in doped semiconductors 
or insulators like Ti02 are the carbon or nitrogen atoms. The magnetic moment is produced 
by the p-electrons which become polarized because of the flat p-bands of these impurity 
atoms. This is called p-electron magnetism and has been investigated intensely during the 
last years. 


II. CRYSTAL STRUCTURE AND DOPING CONFIGURATIONS 

There are three common polymorphs of Titanium dioxide: rutile, brookite and anatase. 
Rutile is the thermodynamically most stable modihcation of Ti02 and is an indirect wide 
band gap semiconductor with an experimental band gap of 3.0 eV.^® 

The rutile structure belongs to the P42/mnm (No. 136) tetragonal space group with 
unit cell parameters of a = 6 = 4.587 A and c = 2.954 A. In rutile Ti02 every titanium 
atom is octahedrally coordinated to six oxygen atoms. The so formed octahedrons show 
an orthorhombic distortion, with the apical Ti-0 bond length being slightly longer than 
the equatorial Ti-0 bond length. Each TiOg octahedron is in contact with 10 neighbour 
octahedrons. The Ti02 rutile crystal structure therefore can be seen as a chain of edge and 
corner-sharing TiOg units. 

Beside the investigation of the electronic properties of pristine Ti02, we examined the 
effect of C and N atoms incorporated in various concentrations and positions into the host 
matrix. In total seven different levels and configurations of doping were examined using 
ab-initio calculations. Fig.l depicts the different positions of the C and N dopants in the 
Ti02 2x2x2 (48 atoms) supercell. In one of the performed calculations we assumed a single 
oxygen atom of the host matrix to be replaced by C or N (position (1)). This corresponds to 
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a doping concentration of 2.1%. Further, two substitutional oxygen sites with two C or two 
N atoms, respectively, were occupied leading to a doping rate of 4.2%. This was done for 
two different conhgurations investigating the effect of increasing distance between the two 
impurity atoms. In the hrst conhguration of this double substitutional state the distance 
between the two dopants was chosen to be minimal (position (1) and (2)). In the second one 
the two dopants were placed on arbitrary oxygen sites (positions (1) and (3)) increasing the 
distance. Moreover, two different interstitial positions of C and N (positions (4) and (5)) in 
the Ti02 host lattice were examined. Finally, we performed a set of calculations exchanging 
these interstitial dopants by their next nearest oxygen neighbour, so that O and C or N 
respectively switch position, with a doping concentration of again 2.1%. 



FIG. 1. (a) and (b) shows the 48-atom supercell of Ti02. The large blue and the small red spheres 
represent the Ti and the O atoms. The orange numbered spheres (1-3) denote the positions of O 
substituted by the C and N dopants and the violet numbered spheres (4-5) notify the interstitial 
doping positions. 


III. CALCULATIONAL DETAILS 

All calculations were performed employing the Vienna Ab-initio Simulation Package 
(VASP)^°^^^ which uses projector augmented wave (PAW) pseudopotentials^® to describe 
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the potential between the ions. The semi-core s and p states for Ti and the 2s states 
for O, C and N respectively were incorporated in the calculations. To save computation 
time, the effects of exchange and correlation were initially treated using the GGA-PBE 
approximation,until all force components were smaller than 0.01 eV/A. During the re¬ 
laxation we allowed for a change of the atomic positions, shape and volume of the cell. The 
hnal convergence was done with the post-DFT HSE06 functional® which is known to yield 
reliable results for the gap size and the position of the impurity bands inside the gap. 

Plane waves with an energy up to 530 eV were included in the basis set, in order to avoid 
Pulay stress and other related problems. The Brillouin-Zone integration was performed 
using a 4 X 4 X 4 P centered fc-mesh with Gaussian smearing set to 0.05 eV. The total 
energy was converged better than 1 x 10“® eV for all cases investigated. The limitations 
of LDA and GGA-I-U in predicting the equilibrium lattice constant and the proper band 
gap^® poses problems in the description of the structural and electronic properties of the 
pristine and doped system. Thus, predicting the position of the impurity states, crucial 
for photocatalytic reasons, becomes only possible after extracting the appropriate U from 
experimental data. 


IV. C/N-SUBSTITUTIONS 

To compare and evaluate the structural changes induced by the impurities calculations 
on pure Ti 02 were done initially to obtain structural and electronic data that can be used as 
a reference. The calculated lattice parameters of pure Ti02 rutile ( a = 4.65A ; c = 2.97A 
) are in good agreement with the experimental values ( a = 4.59A ; c = 2.95A ). The 
direct band gaps calculated using the HSE06 functional ( 3.24eE ) show a considerably 
improvement over the GGA results ( 1.77eV and are again in in good agreement with 
experimental values ( 3.1eV 


A. Single Carbon and Nitrogen Substitution 

In total three substitutional cases were investigated, i) one single oxygen becomes substi¬ 
tuted by G or N impurity. Further, two oxygen sites are substituted by either two carbons 
or two nitrogens, considering two different cases ii) minimizing the distance of the impuri- 
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ties occupying two adjacent oxygen sites, and iii) the impurities occupy two distant oxygen 
positions. The latter two cases allow to study the dependence of the magnetic coupling on 
the distance of the magnetic ions. 



energy [eV] 


FIG. 2. Density of states of the non-magnetic (upper panel) and magnetic (lower panel) state of 
C snbstituted Ti02. The DOS shown includes only the NN aronnd the C impurity, thus 2 Ti and 
5 O atoms. 

Substituting just one oxygen atom by a carbon impurity leads to a magnetic moment 
of 2/is per supercell and arises mainly from the C atom with 0.717/iB per carbon atom 
(measured within the muffin tin radius 1.63 A). Substituting oxygen by carbon introduces 
two holes into the system. As a consequence one complete band on the average becomes 
unoccupied leading to a semiconducting ground state but with a magnetic moment. The 
total moment of per unit cell is composed of the magnetic moment at the carbon site 
plus the polarizations of the surrounding oxygen. Figure 2 depicts the density of states 
(DOS) for the single substitutional case. The upper panel shows the non-magnetic case of 
carbon substitution. The lower one depicts the DOS of the equilibrium magnetic state with 

The substitution of N introduces one hole and consequently leads to a magnetic moment 
of l/is per supercell where 0.57/^5 are located at the N site (measured within the muffin 
tin radius 1.40 A). This result agrees with the experimental hnding of room temperature 
ferromagnetism in N-doped rutile Ti02 hlms by Bao et al.^^ who report a nitrogen moment 
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of about 0.9/is. The DOS for the magnetic and the non-magnetic state is shown in Figure 

3. 



FIG. 3. Density of states of the non-magnetic (upper panel) and magnetic (lower panel) state of 
N substituted Ti02. The DOS shown includes only the next neighbours around the N impurity, 
thus 2 Ti and 5 O atoms. 

Introducing dopants into the Ti02 rutile crystal leads to a change in the atomic distances. 
Changes of the bond lengths between the oxygen (of the undoped system), the substituted 
carbon and the substituted nitrogen atom to its next neighbours (NN), respectively, were 
compared. The atomic distances of the pristine crystal lattice were normalized to those 
of the doped systems using its relaxed volume. By doing so we considered increased dis¬ 
tances caused by the larger ionic volume of the C and N dopant. Table I gives an overview 
on the bond lengths in the respective system. Both the N and the C atom causes an en¬ 
largement of the atomic distances in the crystal. Due to the higher electronegativity of 
N(3.04)>C(2.55),^^ nitrogen tends to interact more with the host semiconductor in partic¬ 
ular with the neighbouring oxygens. This leads to shorter bond lengths for N to its next 
neighbours than for C. These results are again in line with the DOS shown in hgure 3 and 
4. The larger bond lengths and the consequently weaker bonding to the host semiconductor 
makes the 2p states of the C atom more localized in the gap and the C states which occur in 
the VB and CB due to the weaker interaction with O are smaller than in the N case. The 2p 
states of nitrogen are mostly located in the VB and CB, induced by the stronger bonding. 
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As a consequence carbon doped Ti02 rutile has a smaller band gap than the nitrogen doped 
system (table I). 

To investigate the stability of the doped system we calculated the defect formation ener¬ 
gies according to the following formula: 

Eform^X] = E,4X] - E,4bulk] + nifio - hx) (1) 

where Etot is the total energy of the doped supercell with one oxygen atom replaced by the 
impurity X(X=C,N), Etot[bulk] the total energy for the pristine system, n indicates the 
number of oxygen atoms that have been replaced in the supercell by dopant atoms, /io and 
fix are the corresponding chemical potentials. 

The stability of the various conhgurations differ with the oxygen chemical potential, fio 
describes the oxygen environment during synthesis and therefore effects the defect formation 
energies. We dehne the oxygen chemical potential being /io=^h 02 +/ao' with fio' ranging 
from 0 eV to -4 eV, whereas /io'=0 eV dehnes the oxygen rich and fio>=-^ eV the oxygen poor 
case. The oxygen poor case is approximately half the HSE calculated formation enthalpy of 
rutile Ti02 with AHf{Ti02) = —9.95eV which is in good agreement with the experimental 
value of 9.80 eV.^^ To give a more conceptual measure of the oxygen concentration we 
convert fio' to oxygen pressure at a temperature of lOOOK (top x-axis), typical for annealing 
of rutile Ti 02 .^® 

In the tables below we list the defect formation energy for the different doping conhgura¬ 
tions considering the oxygen rich case. The chemical potentials were calculated with respect 
to C in Diamond, N in N 2 and O in O 2 . In addition phase diagrams and E-f°^^[X]/ fio> plots 
are presented in Fig.9. 

B. Multiple Carbon and Nitrogen substitution 

As a second step we study the substitution of 2 oxygens by 2 carbons or 2 nitrogens at 
different but crystallographically equivalent sites. The results are given in table II and the 
positions of the substituted atoms correspond with the site numbering given in Fig.l. For all 
4 cases investigated we hnd a magnetically ordered ground state with Afis per supercell for 2 
carbons and 2/is per supercell for 2 nitrogens. The additional substitution leads to a further 
increase of the cell volume, but considerably smaller than for the single impurity (compare 



TABLE I. Single carbon and nitrogen substitution: Cell volume, bond lengths of ^-"Ti^guatoriai/axial 
and X-0[2], total magnetic Moment Mtot of the supercell, defect formation energy energy 

difference of the non spin polarized and spin polarized state /S.E[NSP — 5P] and band gap for 
spin-up and spin-down. 


x= 

C 

N 

0 

Cell Volume (A^) 

519.95 (-M.9%) 516.18 (-fl.1%) 510.06 

X-Xigqiiatoriai (A) 

2.07 (-h2.8%) 

2.02 {+0A%) 

2.00 

X-Tiaxiai (A) 

2.13 (-h7.9%) 

2.02 {+2.Q%) 

1.96 

x-0[2] (A) 

2.74 (-F6.0%) 

2.62 (-M.9%) 

2.57 

Mtot (iJ'b) 

2 

1 

0 

Eform^X] (eV) 

9.42 

5.43 

- 

AE[NSP - SP] (eV) 

0.24 

0.80 

- 

Spin up gap (eV) 

2.20 

3.20 

3.24 

Spin down gap (eV) 

2.06 

2.60 

3.24 


to table I). The formation energy per atom remains almost constant and in general shows a 
slight increase with respect to the single impurity, only in the case of 2 distant nitrogens a 
very small reduction is found. Studying 2 impurities gives us the opportunity to discriminate 
between a ferromagnetic (FM) and an anti-ferromagnetic (AFM) coupling. In general we 
find the AFM state to be lower in energy, only for the case of 2 distant nitrogens an extremely 
small FM stabilization energy has been calculated, however, this small energy change is at 
the verge of the numerical accuracy. As expected, we hnd that for neighbouring impurities 
(C(l,2) and N(l,2)) the coupling energy is at least one order of magnitude larger than for 
the distant ones (C(l,3) and N(l,3)) which appear to be essentially decoupled. The largest 
energy gain is found for the case C(l,2) where our calculation shows, that in the FM case the 
spin-down gap almost vanishes, so that the opening up of the gap in the AFM state leads 
to the observed energy gain. The band gaps are reduced compared to the case of a single 
impurity, the strongest effect is found for C(l,2) where the interaction of the neighbouring 
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C atoms reduces the gap to 1.4 eV. For the N substitution we hud that the ’’distant” N(l,3) 
case is lower in energy, while for the C substitution the ’’close” C(l,2) conhguration is more 
stable, which is again a consequence of the stronger interaction and points to a tendency of 
the C atoms to form clusters. The multiple nitrogen substitution conhguration N(l,3) has 
a lower per impurity atom than the single substitutional case N(l). 

TABLE II. Multiple carbon and nitrogen substitution: Cell volume, total magnetic Moment Mtot 
of the supercell, defect formation energy per impurity atom energy difference of the FM 

and AFM state AE[FM — AFM] and direct band gap for spin-up and spin-down. 


x= 

C (1,2) 

N (1,2) 

C (1,3) 

N (1,3) 

Cell Volume (A^) 

522.85(-f2.5%) 519.51(-bl.8%) 522.85(-f2.5%) 517.87(-f 1.5%) 

Mtot (pb) 

4 

2 

4 

2 

Eform^X] (eV) 

9.62 

5.47 

9.65 

5.41 

AF[FM - AFM] (eV) 

0.160 

0.018 

1.14x10-3 

-0.019x10-3 

Spin up gap (eV) 

1.40 

2.45 

1.90 

3.16 

Spin down gap (eV) 

1.40 

2.45 

1.90 

2.33 


V. C/N-INTERSTITIAL 

The hnal investigation deals with carbon and nitrogen placed on interstitial sites (sites 
(4),(5) in Fig.l) and their interaction with the respective neighbouring oxygen. Upon putting 
carbon on either of the 2 interstitial sites we hnd that the electronic structure is very similar 
and no magnetic moment appears. We thus restrict ourselves to present only the density of 
states for interstitial site (4) (see hgure 5), the respective electron density is given in hgure 
6(a). The interaction of the interstitial carbon with the host lattice is almost negligible 
and C forms atomic like hat bands at the bottom of the Ti02 gap, which is also rehected 
in the electron density. The non-magnetic state is easily explained from the even number 
of electrons. In contrast to free carbon, which has a total spin of 1, the small but present 
crystal held from the host lattice leads to a double occupation of a single p-orbital and 
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(a)C(l,2) (b)C(l,3) 




energy [eV] 


(c)N(l,2) 


(d)N(l,3) 


FIG. 4. Densities of states for the substitution of 2 oxygens by 2 carbons or 2 nitrogens. The 
densities of states are given for the calculated ground state which is AFM for C(l,2), C(l,3), 
N(l,2), and FM for N(l,3). The ’’bulk” oxygen and titanium states are plotted for 5 neighbouring 
oxygens and 2 titaniums. 


hence a total spin of zero. The same mechanism can be applied to nitrogen with the only 
difference that N has 3 p-electrons and consequently always one singly occupied orbital with 
a resulting total spin of Energetically we hnd that C/N sits on a saddle point of the total 
energy surface. Upon shifting C/N slightly to either side, the interstitial atoms approach 
their oxygen neighbours and form CO and NO molecule type entities which is accompanied 
by a dramatic change in the electronic structure. The highly reactive oxygen tries to reach 
its 2“ state, by forming a CO or NO molecule which ends up almost at the original oxygen 
position. The electron density once the C(5)-0 dimer formed is shown in hgure 6(b). We 
also performed a second set of calculations, where C/N are placed on an oxygen position and 


11 



























energy [eV] 


FIG. 5. Densities of states for interstitial carbons or nitrogen placed on site (4). The densities of 
states are given for the calculated ground state which is non-magnetic for carbon, and magnetic 
for nitrogen. The ’’bulk” oxygen and titanium states are plotted for 5 neighbouring oxygens and 2 
titaniums. 

this oxygen is put at an interstitial site. Starting from this conhguration and after relaxation 
again led to the formation of CO and NO dimers with the same total energy as before. The 
only difference is that the C/N interstitial sits at an energy saddle point, while the oxygen 
interstitial immediately starts to move towards the C/N. Fig.7(a) and (c) shows the charge 
densities for the CO and NO case. One clearly sees the molecular structure with a bond 
length of 1.22 A for CO and 1.31 A for NO. For carbon the bond length agrees with the 
CO double-bond length of 1.20 A, for nitrogen the agreement is less good and lies between 
the NO single-bond length of 1.45 A and the double-bond length of 1.17 A^^ which may be 
caused by the fact that the NO ’’molecule” exhibits a magnetic moment of l/is. 

Table III summarizes the results of our calculations. The hrst 2 lines contain the data 
for the C/N interstitials at their saddle point position. Also for this non-equilibrium state a 
semiconductor is found. Again nitrogen becomes magnetic with l/is, while carbon remains 
non-magnetic. The stable conhguration is found, when both carbon and nitrogen become 
inserted on site (5). For both cases a CO or NO dimer is found, which considerably lowers 
the total energy. In an experimental study done by Chen et al.^® nitrogen doped Ti02 
nanocolloids have been prepared and their photocatalytic properties were investigated. XPS 
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measurements done by this group suggest the formation of NO binding regions in these 
samples, which again underlines our hndings of the NO dimers. Asahi et al.^® report about 
experiments supported by ab-initio calculations where they study various kinds of N complex 
species incorporated into a Ti 02 anatase host matrix and as well encountered the formation 
of NO bonds. DiValentin et al.^^ again found single bonded CO dimers from ab-initio 
calculations. From our calculations we suggest the formation of NO and CO dimers in 
Ti 02 , a possible experimental proof could be the detection of the (symmetric) stretching 
mode frequency, which we calculate to be 1605cm“^ for CO and 1234cm“^ for NO inside 
the Ti 02 lattice. 

The defect formation energies are lower for the interstitial conhguration compared to 
the oxygen substitutional case. For N(5) decreases 1.5 eV compared to the single 

substitutional case N(l) and for interstitial C(5) and substitutional C(l) even differ 

by 4.35 eV. 

TABLE III. C/N-interstitial: Atomic distance of C/N and its nearest neighbouring oxygen, total 
magnetic Moment Mtot of the supercell, defect formation energy per impurity atom and 

direct band gap for spin-up and spin-down. 


C(4) N(4) C(4)-0 C(5)-0 N(4)-0 N(5)-0 


Atomic Distance C/N-0 (A) 1.80(+6.4%) 1.81(+6.9%) 

1.22 

1.24 

1.32 

1.32 

Mtot {fJ-s) 

0 

1 

2 

0 

1 

1 

(eV) 

6.92 

5.33 

6.04 

5.05 

4.08 

3.91 

Spin up gap (eV) 

2.59 

3.51 

0.77 

1.92 

2.89 

3.16 

Spin down gap (eV) 

2.59 

3.40 

3.39 

1.92 

2.32 

2.35 


The densities of states (Fig.8) show the basically different behaviour as compared to the 
C/N impurity on the saddle point state (see Fig.5). The oxygen states are below the bulk 
oxygens and the C/N p-states are inside the gap. The charge transfer of 2 electrons leads 
to a strong splitting of the occupied oxygen and the unoccupied carbon/nitrogen states. It 
is noteworthy that there exists also an interaction of the high lying carbon p-states and the 
Ti d-states. The NO interaction is weaker, which is established by the larger bond length. 
The oxygen states are again below the oxygen bulk, but the nitrogen p-states are just at the 
bottom of the Ti02 gap. The unpaired electron leads to a magnetic moment of Ips and the 
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(a)C(5) (b)C(5)-0 


FIG. 6. Electron-density on a logarithmic scale for substitutional C and the CO dimer in the Ti02 
lattice, lighter colors denote larger electron densities, (a) Electron-density of the non-magnetic 
saddle point configuration of C; (b) Electron-density of the stable configuration of CO. (b) and 
(a) are different planes of the crystal lattice, since the bond-axes of the CO dimer turns during 
formation out of plane (a). 

interaction with the Ti d-states is much smaller than for the carbon case. The respective 
charge and spin density is shown in Fig.7(c), (d). 

For the case of carbon with oxygen on position (4) we find a metastable magnetic solution 
with a magnetic moment of 2fiB per supercell. Unlike the cases discussed earlier in this paper, 
the magnetic moment is no longer located at carbon only, but is distributed between carbon, 
oxygen and 3 neighbouring Ti atoms. Figure 7(a), (b) shows the electron density and the 
spin density, respectively, where the polarization of the Ti d-states is easy to recognize. 
However, due to the unfavourable total energy of this magnetic state we expect that this 
case will not be present in real samples. 


VI. DISCUSSION 

Fig. 9 shows the calculated and ij^n/ l^o' phase diagrams as well as the correspond¬ 
ing defect formation energy versus oxygen chemical potential po'. In Fig.9(a) pc 

vs. no' is plotted giving the thermodynamically stable phases of carbon doped Ti 02 de¬ 
pendent on the chemical potentials of C and O. On the top x-axis the oxygen chemical 
potential is converted into equivalent oxygen pressure at a hxed temperature (T=1000K). 
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(c)N(4)-0 (d)N(4)-0 


FIG. 7. Electron (logarithmic scale) and spin-density (linear scale) for substitutional CO and NO 
dimers in the Ti02 lattice, lighter colors denote larger electron densities. Electron (a) and spin 
density (b) of the metastable magnetic (2/rs) configuration of CO; Electron (c) and spin density 
(d) of the stable magnetic (l/^s) configuration of NO. 

Highlighted in color are the thermodynamically different phases for the single carbon sub¬ 
stitutional case C(l). The phase boundaries for the multiple substitutional cases C(l,2) 
and C(l,3) are indicated by grey and black lines, respectively. Comparing these three phase 
boundaries shows that they are almost identical over a large range of po' and only differ 
for very small chemical potentials of oxygen. For a large range of po' and nc the pure 
undoped Ti02 is the thermodynamically preferred phase (Fig.9(a)). Only for very small 
oxygen chemical potentials, i.e. an oxygen poor phase, the substitution of oxygen by carbon 
becomes favourable. Interstitial doping of rutile is energetically practicable only for higher 
carbon chemical potentials together with an oxygen rich environment. Different values of 
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FIG. 8. Densities of states for the interstitial C(5)-0 and N(5)-0 doping configuration. The 
densities of states are given for the non-magnetic state for carbon, and the magnetic state for 
nitrogen. The ’’bulk” oxygen and titanium states are plotted for 5 neighbouring oxygens and 2 
titaniums. O* is the oxygen of the CO or NO dimer. 

correspond to different reservoirs of carbon atoms, pc = corresponds 

to a reservoir of C in Diamond and pc = l^co 2 —/io2=-15.17 eV to C in CO 2 . That indicates 
that for interstitial carbon doping using CO 2 as a C reservoir only the substitutional and 
the pristine crystal phases are stable configurations. Fig. 9(b) shows the intersection points 
of the defect formation energy for the different doping configurations C(l), C(l,2), C(l,3). 
As a reference for the carbon chemical potential Diamond was chosen as a reservoir. The 
grey line is the interstitial carbon doped case C(5) which has a constant defect energy and 
does not vary with the oxygen chemical potential. Only for /ro' about < —4.5eV the sub¬ 
stitutional phases become thermodynamically more stable. Since all intersection points for 
the different doping configurations are located close to each other at extremely low oxygen 
pressure, sole thermodynamical reasoning would only allow for interstitial doping. 

Fig. 9(c) shows the hn over po/ phase diagram. In contrast to Fig. 9(a) the thermo¬ 
dynamically stable are of pure Ti02 is smaller than in the carbon case. Thus, looking at 
the corresponding oxygen pressure and at the nitrogen chemical potential, all three doping 
configurations become feasible. Nitrogen atoms taken out of a NO 2 reservoir correspond to 
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(a)Phase diagram /ic/^O' 


(b)C-doped /fiQ' 


O 2 pressure at 1000K (atm) 



Og pressure at 10OOK (atm) 

10 '^° 10 '’® 10 '^° 10 '^ 1 10 ® 



(c)Phase diagram ^Ar//io' 


(d)N-doped /fio' 


FIG. 9. (color online)(a) and (c): phase diagrams as a function of the chemical potentials 
and /Tat; (b) and (d): formation energies as a function of the oxygen chemical potential 

llO'- 


a chemical potential of /iat = 11 NO 2 ~ /io2=-ll-36 eV and a N 2 reservoir leads to a nitro¬ 
gen chemical potential of hn = |/rAr2="10-249 eV. Again, as in the CO 2 case, NO 2 would 
only allow for a stable substitutional and pure phase at varying oxygen chemical potential. 
Fig.9(d) shows the intersection of the interstitial and substitutional doped defect formation 
energy lines, where N 2 is assumed as the nitrogen reservoir. In contrast to the carbon case, 
for nitrogen doping the intersection point are in a region of realistic oxygen pressures, which 
may allow for both substitutional and interstitial doping. 

An efficient photocatalyst is supposed to have a ’’clean” band gap with no additional 
states inside that can act as electron-hole traps and recombination centers and therefore 
reduce the lifetime of the photo-induced charged carries, leading to a decrease of photocat- 
alytic efficiency. By doping our system with C and N atoms we try to narrow the band gap 
in an efficient way without introducing mid-band gap states which cause the recombination 
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rate of the charge carriers to become maximum, as derived by Shockley, Read and Hall.^^’"^^ 
Interpreting our calculated DOS and keeping the SRH model in mind, we see that even if 
there is a reduction of the band gap there can be decreased photocatalytic efficiency. 

Doping the Ti02 rutile bulk structure with nitrogen atoms does not lead to the creation of 
impurity states at the mid-band energy level {Emb - dehned as half of the energy difference 
of the top energy level of the VB and the bottom energy level of the CB). In the cases 
studied here, the dopant levels created inside the band gap are mostly located in the lower 
or top third of the band gap. Hence the life time of the charge carriers is reduced but does 
not reach its minimum like for the case of impurity states at Emb- Comparing the different 
nitrogen doped systems, one sees that the spin-up gap remains almost unchanged, showing 
a maximum reduction of 0.08 eV. The spin-down gap on the other hand shows a minimal 
reduction of up to 0.91 eV. For the case of single nitrogen doping there is a negligible small 
narrowing of the spin-up gap by 0.04 eV. The spin-down gap is reduced by a dopant state in 
the upper third of the gap by 0.64 eV. For the case of two nitrogen impurities the spin-down 
gap again is rednced by 0.79 eV for the impurity atoms located at (1) and (2) and by 0.91 eV 
at (1) and (3). For the first case the spin-np gap again is reduced by 0.79 eV and for the latter 
one narrowed by 0.08 eV. In both systems there is, like in the case of single nitrogen doping, 
an impurity level introduced in the upper third of the band gap near the CB. Regarding the 
unstable interstitial state we see a blueshift since both gaps are widened. The spin-up gap 
by 0.27 eV and the spin-down gap by 0.16 eV. For the case of the NO bond like formation 
the band gap narrowing for the spin-up state again is inconsiderably with a rednction of 
0.08 eV. The spin-down gap is reduced by 0.89 eV. Impurity states are introduced in the 
npper and lower third of the band gap. In total nitrogen doped Ti02 causes a redshift in 
the absorption and introdnces band gap states that do not cause the recombination rate to 
maximize and therefore conserve a reasonable photocatalytic efficiency. The only case where 
a blneshift could be detected was in the unstable case of interstitial doping for nitrogen at 
the positions (4) and (5) without the formation of NO dimers. Hence, there have to be other 
reasons leading to the blueshift observations in N-doped rutile crystals. 

Doping with carbon always leads to a signihcant redshift but introduces mid-band gap 
levels at Emb, thus maximizing the recombination rate of the electron-hole pairs causing 
minimal photocatalytic efficiency. For the case of single carbon doping the spin-up band 
rednces by 1.04 eV and the spin-down gap by 1.18 eV. New states are created inside the 
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band gap but not in the region of Emb- For the case of two carbon dopants at the positions 
(1) and (2) mid-band gaps are created and the gaps for both spin directions are reduced 
to 1.40 eV which is beyond the visible spectrum, rather unfavourably for photocatalysis. 
Carbon at the substitutional positions (1) and (3) in contrast shows a desirable behaviour 
with no impurity states in the region of Emb and a band gap of 1.90 for both spin states. 
Carbon on the interstitial position (4) and (5) is similar to the latter case, there are no states 
in the mid-band gap and a gap reduction to 2.59 eV. In the case of the CO dimer there 
are again states at the energy level Embi making the system unattractive for photocatalytic 
application even with a putative band gap of 1.92 eV. In summary there are only two 
favourable conhgurations of carbon in the rutile Ti 02 lattice, the single carbon impurity 
and the one with the carbon impurities at positions (1) and (3). When we replace oxygen 
by C/N both impurities become magnetic. For two impurities on oxygen sites we hnd an 
anti parallel orientation of the magnetic moments to be favoured, however, the low energy 
gain due to the formation of magnetic order suggests that the ordering temperatures will 
be rather low. When C/N is placed on an interstitial site, C remains non-magnetic while N 
always shows a magnetic moment of l/ie- Once the CO and NO dimers are formed, nitrogen 
remains with I^ib while carbon is non-magnetic. Test.^^^^® 


VII. CONCLUSION 

The magnetic structure of carbon and nitrogen doped rutile Ti02 was investigated for 
different doping conhgurations and concentrations. Carbon shows to induce a magnetic 
moment of 2/iB per unit cell for the substitutional doping case and a nonmagnetic solution 
for the interstitial doping conhguration where a CO dimer is formed. For both conhgura¬ 
tions the bandgap is reduced, whereas for the CO dimer case midband states are created 
which could lead to a reduction in the photocatalytic performance. Two substitutional C 
impurities close to each other show AFM coupling and for the case of maximum distance 
paramagnetic behaviour. Relating to the phase diagram carbon is most stable as an inter¬ 
stitial impurity for low chemical potentials. Nitrogen induces a magnetic moment of l/ie for 
both the substitutional and the interstitial doping conhguration. For latter a NO dimer is 
formed like in the interstitial Carbon doping case. Nitrogen doping reduces the bandgap of 
rutile Ti02 as well, whereas no midbandgap state is induced as for the Carbon interstitial 
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doping case. N doping should therefore lead to an increase in photocatalytic efficiency which 
again has shown by several experiments.Looking at the calculated phase diagram, there 
is also a wide range of partial oxygen pressures where substitutional N doping shows to be 
energetically preferable, in contrast to Carbon. The magnetic coupling of two adjacent N 
atoms show to be AFM and like in the Carbon case and for to N atoms with greater dis¬ 
tance to be paramagnetic. Paramagnetic behaviour for the N doping case, however, has 
already been observed in experiments^ and can be reproduced by our calculations. Recent 
experimental groups^^’S^’"^® which investigated N-doped rutile single crystals, indeed mea¬ 
sured ferromagnetic behaviour but cannot exclude the role of vacancies , created by the N 
ion implantation technique as possible sources of ferromagnetism. 
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